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I. INTRODUCTION 



In almost all modern cosmological models, galaxies, clusters, and all large-scale struc- 
tures develop through the gravitational instability of small-amplitude, seed density fluc- 
tuations. In most of these models cold dark matter is an important constituent of the 
total mass density of the Universe. There are two basic types of seed density fluctu- 
ations, curvature and isocurvature|] and in general, both are expected to be produced 
in the early Universe. By definition, the total energy density in an isocurvature fluc- 
tuation is constant; the fluctuation is in the relative contribution to the total energy 
density of different components in a multicomponent system. Important examples of 
this type are the fluctuations induced in the baryons by some dissipative process in a 
Universe containing both baryons and dark matter, and topological or non-topological 
field configurations like cosmic strings or textures. While the amplitude of either type 
of fluctuation on large scales is strongly restricted by microwave background anisotropy 
constraints, the amplitude of small-scale fluctuations can be large, even non-linear, at 
the epoch of last scattering. The spectrum of small-scale fluctuations do not necessarily 
have to reflect the shape of the power spectrum of the primordial fluctuations generated 
at the inflationary epoch, since the small-scale fluctuations may well be generated later, 
e.g., during various cosmological phase transitions. 

In this paper we are interested in isocurvature fluctuations that enter the horizon be- 
fore the temperature of equal energy densities of matter and radiation, Teq = 5.5 VL^h? eV 
0. We will consider scales much smaller than the horizon, so the radiation energy density 
should be homogeneous. 

It is well known || that the growth of small-amplitude isothermal fluctuations is 

1 The division of fluctuations into curvature and isocurvature is strictly true only outside the Hubble 
radius 0. 
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suppressed by the cosmological expansion, and the fluctuations do not grow until after 
the equality epoch. However, this is true only in linear theory. The self gravity of large- 
amplitude, non-linear fluctuations may become important before T E q, and consequently 
they collapse earlier. Therefore they are capable of producing very dense objects after 
they separate out from the general expansion and virialize. 

We refer to these isothermal fluctuations as "clumps. "0 Let us specify the density of 
a dark-matter clump as 

Vdm/pdm = (1-1) 

where $ is not necessarily small. For example, "typical" axion miniclusters Q will have 
$ ~ 1. In Ref. || it was found that accounting for non-linear effects in the evolution of 
axions at the crucial epoch when the axion mass switches on can lead to considerably 
larger density in many miniclusters, with $ in the range 1 to 10 4 . Dark-matter clumps 
seeded by wakes induced by cosmic strings or by textures also will have $ ~ 1 || . Seeded 
clumps are particularly interesting in the case of WIMP dark matter 

It was pointed out in Ref. that the final virialized density in a clump has to scale as 
p F ~ $ 4 p E q. Because of the dependence upon the fourth power of $, even a small increase 
in $ is very important. For the same reason, the final density can be sensitive to the 
details of the evolution of the clump in the radiation-dominated era. To our knowledge, 
a detailed study of the non-linear evolution of large-amplitude isothermal fluctuations 
has never been performed. However, it is very important in various phenomenological 
implications including both direct and non-direct dark matter searches. In this paper we 
consider this problem. 

The dumpiness of the dark matter has important implications for attempts to detect 

dark matter. Clearly the signal in direct detection experiments for dark matter is pro- 
vince the clumps may be very dense compared to the background, we do not refer to them as 
"perturbations." 
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portional to the dark matter density. For the rare direct encounter with a clump, there 
could be a huge amplification of the signal. However, if the dumpiness is too high, the 
flux of unclumped dark matter will be too small for a reasonable detection rate. The rate 
of WIMP annihilation contributing to the 7-ray background 0|J is proportional to the 
density as well. In the case of clumped dark matter, there will be stronger constraints 
on dark matter from indirect searches. In very dense axion clumps, Bose star formation 
becomes possible || (clumps with $ ~ 30 already satisfy the critical condition for this 



||), which in turn can lead to the formation of radio sources []TpJ. Another possible 
manifestation of high density clumps is the phenomenon of microlensing. 

To study the structure and evolution of high density clumps, a full 3-dimensional 
numerical simulation is needed. However, for an isolated clump some relevant physical 
information can be extracted from a 1-dimensional spherical model. The spherical model 
proved useful in studies of the gravitational non-linear evolution in the epoch of matter 
domination when it is possible to find exact analytic solutions [TT| . In the present paper 
we generalize this model to include radiation. Although there are no analytic solutions, 
the result turns out to be very simple: The final density in a virialized clump is pf ~ 
140$ 3 ($ + 1)peq in the whole range of possible values of $, both for $ ^ 1 and for 
$< 1. 



II. A SPHERICAL MODEL 

Let us consider a spherical region of radius r containing an overdensity of pressureless 
matter in an expanding Universe. In a spatially flat Universe, every overdense region 
eventually reaches some maximum size and re-collapses. The total mass of matter in the 
region inside r, Mtot, is an integral of the motion so long as the region expands. Since we 
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will consider scales much smaller than the Hubble radius, we can consider the radiation 
to be homogeneous, with its time evolution determined by the general expansion of the 
Universe, and not by the local conditions. 

The equation of motion for the radius of the region is 

r = ^PR — ■ (2- 1 ) 

It is convenient to change to the conformal time coordinate, drj = dt/a(t), and then to 

rewrite this equation of motion in the comoving reference frame, r = a(rj)R^(rj)^, where 

£ is the comoving label of a given shell and R^{rj) measures the deviation of the shell 

motion from the uniform Hubble flow of the background Friedmann Universe. In what 

follows we shall omit the subscript £ on R(rj), but it should be understood that there is 

a separate evolution for each shell. 

We shall assume that the scale factor a(r]) satisfies the Einstein equations for an 

Qq — 1 Universe filled with radiation and pressureless matter: 

, 2 8nG 4 .. 4nG o 

a = —^-(pm + Pn)a ; a = —^-p M a , (2.2) 

where prime denotes d/drj. We parametrize the radiation and matter energy densities as, 
Pr = Peq(oeq/o) 4 and pu = Peq(oeq/o) 3 - The solution to the background equations, 
Eqs. Q, is 

a(r ] )=a EQ [2(r ] / V *) + (r ] / V *) 2 }, (2.3) 

where i]~ 2 = 27rG , p E QaEQ 2 /3. 

The equation of motion [Eq. ( |2.1| )] in terms of conformal time is 

aR" + a'R' + - ^aW*) = • (2.4) 

The radiation energy density does not enter this equation explicitly, but its effect is 
encoded in the evolution of the scale factor. We also parametrize the total mass of 



matter inside the shell in terms of the excess over the homogeneous background, denoted 
as $(£) = Spm/pm- The total mass within the region is 

4:71 

M TOT = yPEQa E Q 3 [l + 3 • (2-5) 

Changing from rj to x = a/a-^Q as the independent variable, we finally obtain 

s d?R ( 3 \ dR 1 (I + $ \ , , 

This equation reduces to the Meszaros equation || in the limit of small deviation of 
the shell motion from the general cosmo logical expansion, R = 1 — 8 and 5 <C 1, if we 
assume no excess in total mass of the matter, i.e., $ = 0, 

x{1+x) is- + { 1 + r)d- x -2 s=0 - (27) 

The latter is hypergeometric equation, and its growing mode is 5 = <5q(1+3:e/2), implying 
the well known result that the growth of small fluctuations is significant only after the 
Universe becomes matter dominated. 

We have solved Eq. numerically, assuming R(x ) = 1 at some early time, Xq 1. 
Note that at small x, the second derivative in Eq. ( |2.6|) can be neglected, and the solution 
with initial conditions fixed at xo = is R = (1 — 3$a;/2) 1/ ' 3 ~ 1 — $s/2, where the 
expansion is justified since the solution is valid only at small x. Actually, this is the 
separatrix, i.e., independently of the initial value of R', all solutions tend to it (provided 
Xq 1). The results of numerical integration of Eq. (|2.6|) proved to be insensitive to 
R'(xq) already at xqQ < 10 -3 . For several different values of $ they are shown in Fig. 1. 

It is possible to find an analytic approximation to R(x) as a power series in x to any 
given order. To third order it is 

$x $ 2 x 2 (8$ 3 - $ 2 )x 3 

R = 1 . (2.8) 

2 8 144 ( J 
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The first three terms in this decomposition provide a good practical fit to the solution. 
This fit is shown in Fig. 1 by the dotted line. The last term in Eq. ( |2.8| ) shows that the 
solution is not simply a function of the product Qx. 

Our main goal is to find the parameters of the fluctuation, i.e., its radius and density, 
at the moment when the fluctuation turns around. For later times the assumption of 
spherical symmetry breaks down; however, we can assume that the radius of a virialized 
gravitationally bound object will be one-half of the turnaround radius, and the density 



inside the object will be eight times larger than the density at turnaround [11]. The 
turnaround time, defined by r = 0, in coordinates of Eq. ( |2.6|) is the solution to the 
equation R + xdR/dx = 0. Up to order 0(x 2 ) the function R only depends upon the 
product $£, so to this order in x the scale factor at turnaround will be given by xta = 
const/$ and to the same order in x, Rta = const. The matter density of a fluctuation at 
turnaround is /?ta = (l/4:7rr 2 )dM/dr, with r = R^a^taC,, so it is appropriate to represent 
the parameters at turnaround as 

x TA = C;/$; p TA = C > EQ |^(l + $)£ 3 ■ (2.9) 

We expect only a weak dependence of C x and C p upon since higher order corrections 
are small. Results of a numerical integration for C x and C p are shown in Fig. 2 and Fig. 
3 respectively. These figures demonstrate that for practical applications we can consider 
C x and C p to be constants in entire range of possible values of $, both $ ^> 1 and $ < 1. 

Let us compare our results to predictions of a standard spherical model, which is 
valid for the matter-dominated epoch, i.e., at small $. The spherical model predicts 
Pta/pb = 97r 2 /16 for the density contrast at turnaround, where pb is the background 
density |TT ]. Using Fig. 2 we can extract the corresponding values of /?b and then calculate 
the standard spherical model prediction for C p . This is presented in Fig. 3 by the dashed 
line. Both dashed and solid lines coincide in the limit of small $, as they must. 
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III. APPLICATIONS 



With the function <&(£) in the general form we have assumed, the effective radius 
after virialization for each shell of a given label £ will be half the turnaround radius of 



the shell.fj Using Eq. (|2.9|) we obtain for the final density profile in a virialized object 

PF ^uo PEQ ^-^(i+^e, (3.1) 

where we have set C p ~ 17. For the core density, this formula gives pp ~ 140$q($o + 
1)Peq, where $o = *£(0). The numerical value of the density at equality is peq ~ 
3 x 10~ 16 (Qh 2 ) A g cm -3 . Now let us turn to a few specific examples. 

A. Axion miniclusters 

Fluctuations in the density of axions can be very high, possibly spanning the range 
1 < $ < 10 4 |6,12]. Even with $ as small as 1, the density in miniclusters which form 



out of these fluctuations can be as much as 10 10 times larger than the local galactic halo 
density g]. Using Eq. (|QD, we obtain p F ~ 9 x l(T 14 fWg cm" 3 for $ = 1. With 
$ 3> 1 this result must be multiplied by <E> 4 /2. 

The typical mass of an axion minicluster corresponds to the total mass in axions 
within the horizon at T ~ 1 GeV when the inverse axion mass is equal to the Hubble 
length: Muc ~ 10~ 9 M Q . The present probability of a direct encounter with a minicluster 
is small, the encounter rate is 1 per 30 million years with $ = 1. Although the signal 



in an axion detector |T3[ from a close encounter would be enormous, it might be a long 
wait with a weak signal between encounters if a major fraction of the axions are part of 
miniclusters. 

There should be some miniclusters with $ in the range 10~ 3 <<&<!. These collapse 



3 Strictly speaking this is not always true, but for the practical applications we will consider below 
the approximation should be reasonable. 
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during the matter-dominated epoch and have a larger radii than those with $ > 1 
which collapsed in the radiation-dominared epoch, so the probability of an encounter 
with a clump with 1 can be larger. From the point of view of direct searches, even 
miniclusters with density contrast of order two times the average with respect to the 
galactic halo density are important. Such miniclusters form just prior to the moment 
of galaxy formation and started with $ ~ 10~ 3 . For $ < 1 the expected time between 
encounters is given in terms of the number density of clumps, n, the geometric cross 
section of the clump o ~ -Rclump^ an d the virial velocity v as 

t= J_^]_PfRclvmp^ (32) 
nav $ pn v 

where pu is the halo density, and -Rclump/^ is the time the Earth spends inside the 
minicluster. The factor of $~ 1 appears because the number density of miniclusters with 
$ <C 1 is suppressed in our model. A minicluster with $ < 1 would require the initial 
misalignment angle 9 (which is uniformally distributed in the range to 2tt) to be finely 
tuned to the mean misalignment angle to an accuracy 56/8 ~ $/2. Using Eq. fl2.9p , we 
finally obtain 

r = $r$ =1 = $ • 3 x 10 7 yr- (3.3) 

Note that the miniclusters discussed so far appear if the axion field is uncorrelated 
on scales larger than the Hubble radius at T ~ 1 GeV. However miniclusters with $ <C 
1 can appear from primordial density fluctuations generated by inflation without the 
suppression factor of If this is the case, then r ~ $ 2 • 3 x 10 7 yr. Since $ would 
be small, this would give a reasonable encounter rate, and the question of formation and 
survival of small-scale clumps within the galaxy is worth further study. 

Another astrophysical outcome of very dense axion clumps can be the possibility of 
"Bose star" formation in axion miniclusters. The Bose-Einstein relaxation time in the 
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minicluster due to axion self-interaction is smaller than the present age of the Universe 
with $ > 30 0,1 . 



B. Accretion by a point mass 

The density profile in the halo accreted by a previously formed clump can be calcu- 
lated in the approximation of secondary infall onto an excess point mass of mass m. In 
this case = m/M, where M = 47tpeqcieq 3 £ 3 is the mass of the background dark 
matter within the shell with the label £. Substituting this into Eq. (|3.1| ) we find 



p F « 140p EQ (m/M) 3 . (3.4) 

This can be translated into pj? as a function of r since M has to be understood as the 
mass of dark matter residing within r. The result is pp oc r _9//4 , the same power law one 



usually obtains for secondary infall in the matter-dominated era [ 14 1 . However, Eq. (|3] 
is valid regardless of the time when collapse actually occurs. 

C. Cosmic strings 

We can apply Eq. (|3.4| ) for clumps of dark matter seeded by loops of cosmic strings 
so long as the peculiar velocity of a loop is sufficiently small. Since the string loop is 
not a point object, this formula breaks down in the region of small M. Namely, when 
the given shell turns around at xta ~ 0.7 M/m, the loop size Is has to be smaller than 
the physical radius of the shell, tta = ^ta-Rta^eq, for Eq. ( |3.4|) to be valid. This 
gives the restriction M > 3m 3 / 2 p _3//4 pEQ 1//4 = Mc, where m = pis- We can consider 
Mc as the mass of the core region. The corresponding maximum density which can 



be achieved in the core is pc ~ 15^/ p 3 peq/M C - This value of the core density could 
be many orders of magnitude larger than the density at equality. However, as we see 
from Eq. ( |3.4|) , pp is much greater than peq only in the case when the mass of the 
string loop is larger than the mass of the accreted dark matter. Consequently, with 
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gradual loop decay due to emission of gravitational radiation, the dark matter clump 
will adiabatically expand and diminish in density. This process of clump expansion will 
continue untill m < M. Since in the gravitational field the product rm is an adiabatic 
invariant for each dark matter particle, where r is the effective radius of the orbit, we 
conclude that in any clump for which M < m initially, the present density will have the 
same order of magnitude, p ~ 10 2 /?eq- Dark-matter clumps seeded by wakes induced by 
long segments of moving cosmic strings or by textures also will have $ ~ 1 (see and 
correspondingly the same virialized density, p ~ 10 2 peq- While this density is sufficiently 
high to be significant in applications like annihilation of dark matter particles, it is too 
small to cause microlensing, as we show below. 

Cold dark matter accretion onto string loops both in the matter and radiation dom- 



inated era was also considered in Ref. [15], however, only in the linearized limit of Eq. 
(PD, i.e., Eq. (ET 



D. WIMP annihilation 

In the case where the dark-matter particle species is a stable weakly interacting 
massive particle (WIMP) such as a very massive neutrino or a supersymmetric particle 
(photino, higgsino, or scalar neutrino), the WIMPs can annihilate, contributing to the 
7-ray flux. This places severe constraints on the dark-matter density near the center of 
the galaxy |7|||. Clumped dark-matter annihilation is even more efficient, and places a 
very strong limit on the clumpines as a function of the WIMP properties ||. 

For the 7-ray flux on Earth from WIMP annihilation in the clump, we can write 

L< ~ Lr{m 2 x ' (3 ' 5) 
where r Q is the distance from the Earth to the clump (r m 8.5 kpc is the distance to the 

center of the Galaxy), mx is the particle mass, and M is the mass of the clump. Since 

the particles are non-relativistic, both in the clumps and at the epoch of cosmological 
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freeze-out of the WIMPs, the thermal average of the cross-section in the Eq. ( |3.5| ) is 
directly related to the cosmological abundance [ |T6|j : 
4 x 10~ 27 

( av ) ~ — o — — cm 3 sec _1 . (3.6) 

If we consider a large region of (possibly) clumpy dark matter, like the galactic core 
or spheroid, we must sum up the fluxes from each individual clump. As a result, instead 
of M in Eq. (|3.5| ) we have to substitute £M tot , where £ is the mass fraction of all clumps 
to the total mass M tot in the region. For example, using p ~ 10 2 peq we obtain for the 
central spheroid (M tot « 1O 8 M ) 

7 7 ss £Q 3 x h 6 m^ cm- 2 sec _1 , (3.7) 

where m 2 o = rrix/20 GeV. This has to be compared to the observed upper limit to the 
7-ray flux in the direction of the Galactic center J 7 £3 4 x 10 _7 cm _2 sec _1 . 

E. Gravitational microlensing 

Two conditions must be satisfied for the clump to cause gravitational microlensing 
First, the mass of the clump has to be in a range near 0.1M Q . Second, the physical 



radius of the clump has to be smaller than the Einstein ring radius, Re = 2yGMd 
where d is the effective distance to the lens (typically d ~ 20kpc). The second condition 
restricts the density of the minicluster to be p > 10 7 peq/ \/M_i, where M_ x = M/O.1M . 
If the lensing object is a clump of non-interacting cold dark matter, it has to be formed 
from a density fluctuation with $ > 20. 

Dark matter clumps seeded by string loops or textures, which were considered in Ref. 
0, are in the appropriate mass range; however, they have $ ~ 1. Axion miniclusters can 
have $ > 20; however, they are too light. While it is possible to invent models where 
both conditions are met for some of the clumps (one example could be an axion model 
with an extremely small, but non-zero, value for the u-quark mass), it is hardly likely 
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that a substantial amount of the dark matter has evolved into clumps capable of lensing. 
On the other hand, anticipating significant numbers of microlensing events (for the first 
positive reports see Ref. flT8"| ) in the future, it is not excluded that some of them could 



be caused by the clumps in such classes of models (especially if collisional relaxation is 
significant). The corresponding light curve will be different from the MACHO event since 
clumps are extended objects. 



When our paper was almost completed we became aware of the paper Ref. [19| where 



Eq. (|2.6| ) was studied in details; however, again for the case which corresponds to $ <C 1. 
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FIGURE CAPTIONS: 



Fig. 1: Numerical solutions to Eq. ( |2.6| ) for several different values of $. The second- 
order fit is shown by the dotted line. 



Fig. 2: The coefficient C x in Eq. (|2.9| ) function of $. 



Fig. 3: The coefficient C p in Eq. (|2.9|) as a function of $. The dashed line is the 
prediction of the standard spherical model (matter without radiation). 
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